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ABSTRACT 

We present results of a new feedback scheme implemented in the Munich galaxy 
formation model. The new scheme includes a dynamical treatment of galactic winds 
powered by supernovae explosions and stellar winds in a cosmological context. We 
find that such a scheme is a good alternative to empirically-motivated recipes for 
feedback in galaxy formation. Model results are in good agreement with the observed 
luminosity functions and stellar mass function for galaxies in the local Universe. In 
particular, the new scheme predicts a number density of dwarfs that is lower than in 
previous models. This is a consequence of a new feature of the model, that allows an 
estimate of the amount of mass and metals that haloes can permanently deposit into 
the IGM. This loss of material leads to the suppression of star formation in small haloes 
and therefore to the decrease in the number density of dwarf galaxies. The model is 
able to reproduce the observed mass-stellar metallicity and luminosity-gas metallicity 
relationships. This demonstrates that our scheme provides a significant improvement 
in the treatment of the feedback in dwarf galaxies. Despite these successes, our model 
does not reproduce the observed bimodality in galaxy colours and predicts a larger 
number of bright galaxies than observed. Finally, we investigate the efficiency of metal 
injection in winds and in the intergalactic medium. We find that galaxies that reside 
in haloes with M v -„ < 10 12 /i _1 Mq may deposit most of their metal mass into the 
intergalactic medium, while groups and clusters at z = have lost at most a few 
percent of their metals before the bulk of the halo mass was accreted. 
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1 INTRODUCTION 

Dwarf galaxies are the most abundant galactic objects 
in the universe and play a key role in galaxy forma- 
tion. Because of their shallow gravitational potentials, feed- 
back from star formation is most efficient in these galax- 
ies and is believed to success fully eject a l arge amount 
of their gas and metals (e.g. [Larson 19741: Dekel fc Silkl 



19861: iMac Low fc Ferrarall 19991: Ispringel fc HernquistfeOOa 
Veilleux. Cecil fc Bland-Hawthorn! 120051 and references 
therein). As a consequence, dwarf galaxies may have a fun- 
damental role in polluting the diffuse intergalactic medium 
(IGM, her eafter) with the products of stellar nucleosyn- 
thesis (e.g.lAguirre et al.ll200ll; iTheuns. Mo fc SchavdlioOll ; 



SpriiiR'cl & Hcrnquist 2003; Bcrtonc, Stoehr . 
Qppenheimer fc Davell2006l ). 



White 2005; 



Observations of galaxies and their gas content in the 
local universe indicate that feedback in the form of galactic 
winds plays an important role in regulating star formation 
llGarnettll2002l:lKauffmann et al.|[2003l : iTremonti et~ai1l2004l ; 
ISalzer et al.l 20051 ). As mentioned above, the impact of feed- 
back is believed to be strongest in dwarfs. It is however 
important to note that feedback can be equally effective in 
temporarily reheating the interstellar medi um and quench- 
ing star formation in mo r e massive galaxies dHeckman et al.l 



20001: ICecil et all l200ll: IRupke. Veilleux fc Sanders! |2002|; 



Shaplev et ail 120031 : IRupke fc Veil leux 2005» . A variety of 



observational data suggests that, for each unit of gas mass 
that is converted into stars, an equal or larger amount of in- 
terstell ar gas can be reheated and ejected from the disk of a 
galaxy llMartinlll999l : IPettini et ai1l2002l : IRupke fc Veilleuxl 
2005; lMartmT l2006). Since the ejected mass may be lost to 
the galaxy and deposited into the IGM, this has a dramatic 
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effect on the total amount of stars that can be formed in a 
galaxy during its lifetime. 

By regulating star formation, feedback also alters 
the properties of the entire galaxy population. The 
shape of the luminosity function (i.e. the space den- 
sity of galaxies as a function of their luminosity) is be- 
lieved to reflect the effects of different feedback mech- 
anisms. Different studies have demonstrated that feed- 
back by galactic winds is responsi ble for shaping the 



faint end of the luminosity function (|White fc Rees 



Dekel fc Silklll986l; iKauffmann. White fc Guiderdoni 



1978 



1993 



Somerville fc Primacklll999l ; iBenson et alj 120031 ). More re- 



cent work has focused on feedback from accreting black 
holes in active galactic nuclei (AGN feedback, hereafter) 
as an important physical mechanis m to suppress the cool- 
ing flows in relatively massive haloes (fPi Matteo et al.ll2005l ; 
ICroton et al.ll2006l ; iSower et al.ll2006l ). 

While some of the gas (and metals) ejected by 
winds may fall back onto the galaxy at a later time 
via galactic fountains or gravitational accretion, a frac- 
tion of the ejected material may be deposited perma- 
nently into the IGM. Metals have been detected in the 
lowest density regions of the universe probed by the 
Lycv forest in the spectra of high-redshift quasars, and 
the IGM metallicity has been estimated to be in the 
range 10~ 3 - 10~ 2 Zp ilSongaila fc Cowielll996l;lRauchl | l998 | ; 
Schave et ai1l2000l; ISchave et al.ll2003l ; lAguirre et al-lhooi 



Simcoe. Sargent fc Rauchl 120041 ) Although most of the in- 



tergalactic metals appear to be associated with galaxies 
|Adelberger et alj 120051 ; IStocke et all 120061), a significan t 
fraction does not show such a correlation ( Pieri et al.l l2006). 
Galaxies remain the best candidates to explain the enrich- 
ment of regions at lar ger overdensiti e s, suc h as filaments, 
groups and clusters (lAguirre et al. 1 I2001I; Frii euns et al.l 
20021; |Pe Lucia et al.ll2004l; iBertone. Stoehr fc White! 120051 ; 



Daigne et alj|2006l ; lOppenheimer fc DavetoOrj ). As for the 



enrichment of the lowest density regions, it is still not clear 
which is the contribution from the first populatio ns of stars 
l|Murakami et alj|2005l ; iMatteucci fc Calurall2005l ). 

In this work, we implement a dynamical treatment of 
galactic winds in a galaxy formation model in order to test 
if this leads to any improvement on previous empirically- 
motivated feedback schemes used in semi-analytic models. 
Such a scheme introduces a further level of sophistication 
in the treatment of feedback. To this end, we use an up- 
dated version of the wind model of lBertone. Stoehr fc White] 
ll2005l) and the Munich g alaxy formation model used in 
|Pe Lucia fc BlaizotJ (|2007l ) (DLB07 hereafter). Our wind 
model entirely replaces the supernovae feedback scheme 
used in the semi-analytic model but it does not affect the 
AGN feedback, wh ich does not depend on star formation 
l|Croton et al.ll2006l). We apply our model to the Millennium 
Simulation (|Springel et al.ll2005T ) and we compare our model 
results to those from the latest Millennium galaxy catalogue 
of DLBof]. 

This paper is organised as follows. In Section [2] we de- 
scribe the new features of our feedback scheme and its im- 
plementation in the Millennium simulation. Section [3] shows 



results for the galaxy luminosity functions, the stellar mass 
function, galaxy colours, and morphological fractions at 
z = 0. In Section |3] we discuss the impact of winds on the 
chemical evolution of galaxies. Finally, in Section [5] we in- 
vestigate the role of winds in recycling gas and metals, and 
in enriching the IGM through the cosmic time. We draw our 
conclusions in Section [6] 



2 GALAXY FORMATION 

In this work we use the Millennium simulation^ to trace 
the evolution of dar k matter haloes in a cubic box of 500 
ft _1 Mpc on a side (jSpringel et al.l l2005h . The underlying 
cosmological model is a ACDM cosmology consistent with 
a co mbined analysis of t he 2dFGRS and first year WMAP 
data (|Spergel et al.i r2003). The cosmological parameters are 
n m = 0.25, n b = 0.045, h = 0.73, Q A = 0.75, n = 1, and 
as — 0.9. The Hubble constant is parameterised as Ho — 100 
/i _1 km s" 1 Mpc" 1 . The simulation follows N = 2160 3 dark 
matter particles of mass 8.6 x 10 s /i -1 Mq. Since we require 
dark matter haloes to be formed by at least 20 particles, 
the minimum halo mass is 1.7 x 10 10 /i — 1 Mq, with a corre- 
sponding baryonic mass of about 3 x 10 9 /i _1 Mq. 

As basis for our model, we use the galaxy formation 
model of DLB07. We refer the reader to the original paper 
for more details about the different physical processes ex- 
plicitely modelled and their implementation. We note that 
the model used in DLB07 retains many prescriptions of the 
previous v ersion of the Munich galaxy formation models de- 
scribed in ICroton et al.l (|2006T ). but differs in a number of 
details. DLB07 use a different construction of the dark mat- 
ter merger trees, modified as to effectively track the most 
massive progen itor history. DL B07 adopt a Chabrier initial 
mass function (IChabrieijl2003T ) and an improved treatment 
of dust attenuation which takes into account the attenua- 
tion of young stars within their birth clouds. In this work, 
we use the same set of parameters adopted in DLB07. The 
feedback parameters 6disk, £haio and the reinco rporation pa- 
rameter 7 e j (see Table 1 in lCroton et alj (|2006T ) are not used 
in our model. Our wind scheme also contains free param- 
eters, which we discuss in the following and summarise in 
Table □ 



2.1 A dynamical feedback model 

The feedback model used here is based on the work by 
IBertone. Stoehr fc White] |2005l) . Galactic winds are treated 
as spherical blastwaves, powered by supernova explosions 
and stellar winds, which expand in a cosmological context. 
The physics of galactic winds is described by the equation 
of motion that governs their dynamical evolution. 

The main difference betwe e n th is analysis and 
that of IBertone. Stoehr fc White] (|2005l ) resides in the 
fully self-consistent treatment of feedback. While in 
iBertoneT Stoehr fc White! l|2005h winds are added "on top" 
of the semi-analytic prescriptions for galaxy formation, here 
we integrate the dynamic treatment of winds together with 



1 The catalogue is publicly available at: 
http://www.mpa-garching.mpg.de/millennium/ 
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the existing prescriptions for star formation, chemical evo- 
lution and gas cooling. This results in a scheme in which 
the ejection and recycling of gas and metals is dictated by 
the dynamic evolution o f the winds, rather than by phe- 



19931: ISomervil 
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Primackl 19991: Kauffmann et al.1 1999: 


Springel et al.l 


2001 


; De Lucia et al.1 2004 ICroton et all 



feedback in semi-analyt ic models has recently been pre- 
sented bv lMonacol l|2004 ). 

We describe the dynamical evolution of a wind as a 
two-stage process: a pressure-driven, adiabatic expansion, 
followed by a momentum-driven snowplough. During the 
first stage, a spherical bubble of hot gas blows out the disk 
of a galaxy and a wind is formed. The expansion of the 
bubble is driven by the pressure of the outflowing gas and 
the equation of motion is expressed by the conservation of 
the bubble energy E^: 

dE * 1 »-■».(*) 



-M w v w i 



dAW A n2 
+ 4tvR 2 ■ 



dt 



-p V + U — p 



GM h 



R 



,) - VoPo } (1) 



R and v s are the radius and the velocity of the shock, M w 
and w w (R) the mass outflow rate and the outflow velocity 
of the wind at the shock radius (see Subsec. I2.2|l . p , P , 
u and Wo are the density, the pressure, the internal energy 
and the outward velocity of the surrounding medium and 
Mh the total mass internal to the shock radius. The term 
dAW/df indicates the rate at which gravitational energy is 
transferred from the swep t-up gas to the other components 
|Ostriker fc McKed \wm ■ The expansion of a wind is de- 
termined by the balance between the forces acting on the 
wind: the injection of energy and momentum from SN ac- 
celerates the expansion, while the combined effect of gravity, 
thermal pressure, and ram pressure slows it down and might 
eventually make the wind collapse. 

When radiative losses become dominant and most of the 
energy transferred to the swept-up gas is radiated away, the 
adiabatic stage is terminated and a thin shell of cooled gas 
forms near the outer edge of the bubble. This typically hap- 
pens when the cooling time of the bubble becomes shorter 
than the age of the wind. When the snowplough stage sets 
in, the expansion of the shell is described by the conservation 
of the wind momentum: 



(mv s ) = 



M w (» w (R) — v s ) — ^f h m - 



R? 



AnR [P + p V (Va 



(2) 



where m is the mass of the shell. 

The only formal difference between these eq uations and 
those used in iBertone" Stoehr fc White! |2005) consists in 
the removal of the free parameter e. This means that a wind 
expanding into a halo sweeps up all the hot gas it encoun- 
ters on its path, instead of just a fraction of it. It is therefore 
equivalent to assuming e = 1. In our model winds are pow- 
ered by all galaxies in a halo and effectively blow "from 
haloes" , receiving energy, momentum and mass from all the 
galaxies in the halo (see sec. I2.3[l . When haloes merge and a 
former central galaxy becomes a satellite, its wind is trans- 
ferred to the central galaxy of the merged halo. If the central 
galaxy is already blowing a wind, the two are combined un- 



der the assumption that mass, energy and momentum are 
conserved. The mass and the metals ejected into the IGM 
by the merged galaxies are added to the ejected component 
of the new central galaxy. 

If the momentum of the outflowing material is not able 
to balance the momentum losses by gravity and pressure 
forces, the wind can collapse and create a galactic fountain. 
The dynamics of the collapsing wind can be traced by Eq. ([2| 
until the shell falls back onto the central galaxy. We assume 
that all the mass in the collapsed wind is instantaneously 
reincorporated into the hot component of the halo. In pre- 
vious schemes (|Pe Lucia et al.ll2004l : ICroton et alj|2006h the 
ejected material is reincorporated into the hot component, 
thereafter being again available for cooling, on a time-scale 
that is proporti o nal to the dynamical time-scale of the halo. 
iDe Lucia et~ai1 l|2004 ) argue that various observational re- 
sults can be reproduced only by models in which the ejected 
material is re-incorporated on a time-scale comparable to 
the age of the Universe. 

If a large amount of energy and momentum is injected 
in a bubble, a wind can overcome gravity and pressure forces 
and expand to a large radius. When the wind reaches the 
Hubble flow, it escapes from the halo and its mass and met- 
als are lost for ever from the galaxy. We define the transition 
between the halo environment and the Hubble flow as the 
point where the velocity of the surrounding medium, defined 
as v — RH (z) — « csc (R) /2, is dominated by the Hubble 
flow term and becomes positive. Here, H (z) is the Hubble 
constant at redshift z and v csc (R) is the escape velocity of 
the halo at radius R, calculate d assuming a NFW profile 
|Navarro. Frenk fc Whitelll996l ). A consequence of this as- 
sumption is that the maximum radius to which winds can 
travel is proportional to the virial mass of the halo. Winds 
outflowing from small haloes escape at smaller radii than 
winds expanding out of larger haloes. In the following, we 
will call "mass in IGM" the wind mass that reaches the 
Hubble flow. This material cannot be re-incorporated and 
therefore does not partecipate again in the evolution of the 
galaxy. 



2.2 Mass and metals in winds 

The evolution and exchange of metals between the dif- 
ferent g alactic compone n ts fol lows the same scheme pro- 
posed in IDe Lucia et al.l l)2004h , modified to integrate self- 
consistently the removal of gas and metals from the cold gas. 
When stars form, metals are produced and instantaneously 
returned to the interstellar medium. Subsequently, they can 
be deposited into the hot gas and ejected component by 
feedback processes. 

During the adiabatic stage, the outflowing mass creates 
the over-pressurized bubble of hot gas. The mass of the bub- 
ble Mb is determined by the conservation of mass law: 



dMb 
df 



M w + 4irR 2 p (v s 



(3) 



where the first term on the right hand side represents the 
mass of cold gas reheated by SN explosions and stellar winds 
and the second term defines the amount of hot gas swept up 
by the wind and accreted during the expansion in the halo. 
In a snowplough, most of the mass accumulates in a thin 
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shell of mass m: 

dm 
~dt 



= Mw 



+ 4%R 2 p (« s — v ) ■ 



(4) 



A fraction M w v s /v w (R) of the outflowing ISM does not 
reach the shell and remains in the cavity behind the shock. 
Assuming that the amount of metals exchanged is propor- 
tional to the corresponding amount of mass (i.e. neglecting 
any metal-loading effect), the corresponding equations for 
metals are: 

= M^Z co m + 4irR 2 p (v s — v ) Z hot , (5) 

where Zcold is the metallicity of the cold gas and Zhot the 
metallicity of the hot gas, and: 



dm z 
~df 



= M„ 



(R) 



4n R? p (v a 



Vo ) 



(6) 

A fraction Z co idM w « s /« w (R) of the outflowing metals re- 
mains in the cavity behind the shell. 

Eqs. [3}{6] imply that the amount of mass and met- 
als ejected from the cold gas are M w dt and M v Z co \ddt 
respectively. In our model, we assume that supernova 
ejecta efficiently mix with t he ISM before leavin g the 
disk of the galaxy, as done in iDe Lucia et~ai1 <|2004h . The 
amount of mass and metals removed from the halo hot 
gas and swept up on to a wind are 4nR 2 p a (v B — v )dt and 
A-kR 2 p (v B — v ) Zhotdt respectively. The removal of hot gas 
from the halo is in principl e similar to the imp lementation 
of the "ejected mass" in the lCroton et all l|2006l ) model, but 
here it is justified on the basis of the dynamical evolution of 
the winds and their sweeping up the gas on the path of the 
expansion. 



2.3 Wind initial conditions and model parameters 

Observations of galaxies in the local and high redshift uni- 
verse suggest that the mass loss rate of wi nds is r o ughly 



verse suggest tnat tne mass loss rate ot wi nds is r ouglily 
proportional to t h e star formation rate dMartinl 19991 : 
iPettini et all 120021 : iRupke fc Veilleuxl 120051 ; iMartirJ l200r£ 
A sim ilar behaviour is predicted by the models of Monaco! 
|2004h and lShu, Mo fc Mad (|2005T l. Only a weak dependence 
of the ma ss loss rate on the virial velocity of haloes is inst ead 
observed (|Martinlll999l : lRupke. Veilleux fc Sander J2005I ). In 
our model, we parametrise the mass loss rate M w as a func- 
tion of the star formation rate M* and of the virial velocity 
of the halo v v i r : 



220 km s" 1 



(7) 



In the above equation, M„ represents the sum of the star 
formation rates of all the galaxies in a halo. The power law 
index a w and the constant wind mass loss rate r\ are free 
parameters. According to this choice, all galaxies in a halo 
contribute to power the corresponding wind and to eject 
mass and metals. The dependence of the wind mass loss 
rate on the halo virial velocity «~ ; " w implies that the ra- 
tio M w /M* is highest in the smallest haloes. We note that 
this definition of the mass loss rate i s equi valent to that of 
the reheated mass in IDe Lucia et a l. (2004). In our fiducial 
model we assume r\ = 2 and ct w = 2. 



Model parameter 


Fiducial value 


n 


2 




800 km s" 1 




2 



Table 1. Summary of model parameters. 



At the galactic radius -R ga i, where energy and momen- 
tum are injected into the wind, the initial wind velocity 
u wo is calculated assuming that the energy injection rate 
in winds _E W = M w i4 is constant and independent of halo 
properties. Given the dependence of the mass loss rate on 
iVir in Eq. [7] the wind initial velocity can be expressed as 



220 km s" 



(8) 



where the wind initial velocity V^, is treated as a free pa- 
rameter and we choose as fiducial value Vw = 800 km s _1 . 
Table [T] summarises the values of all the free parameters for 
our fiducial model. The wind velocity at a radius r is cal- 



culated as v w (r) — v. 



2 

wo 



i (r) — 4> (-Rgai)], where <f> is the 



gravitational potential of the halo, calculated assuming an 
NFW profile. 

Observations of galaxies both in the local and 
the high redshift universe estim ate wind velocities 
in the range 100-1500 km s" 1 (IStrickland fc Stevens! 
| 2000l:lFrve. Broadhurst fc Benitej|2002l;IShaplev et alj | 200a 
iRiroke fc Veilleuxl 120051: IRupke. Veilleux fc Sanders! booil : 
see also Veilleux. Cecil fc Bland-Hawthorn| [2005 for a com- 
plete review of the phenomenology of galactic winds) . Since 
winds are multiphase phenomena, the detection of the out- 
flowing gas in different wavebands gives information about 
wind gas components at different temperatures. These may 
yield very diverse values for the wind velocity. The X-ray 
emitting gas is believed to represent the most energetic 
phase of the wind, containing most of the metal-r i ch su - 
pernova ejecta. Simulations bv lStrickland fc Stevens! (|2000t ) 
demonstrate that this gas is outflowing with the highest 
velocities, but it is not representative of the bulk of the 
mass in winds. The latter - that is what our wind model 
is meant to describe - can be detected in less energetic 
wavebands and has much lower veloci t ies (|Rupke fc Veilleuxl 
120051 : iHeckman et alj|2000h . iMonacol (|2004 ) implemented a 
feedback scheme in which the cold and the hot outflowing 
phases of a wi nd are followed individually. This sc heme has 
been used by iMonaco. Fontanot fc Taffonil [|2007r i in their 
galaxy formation model MORGANA. 

Although the large wind velocities given by Eq. [8]seem- 
ing to indicate that winds can escape from large haloes, this 
never occurs in practise. As we will show in the following, in 
haloes with M v i r > 10 12 /i _1 Mq the injected energy is never 
large enough to overcome the gravitational potential and the 
pressure forces of the intra-cluster gas. In these haloes, the 
wind usually collapses within a few time-steps and its mass 
and metals are reincorporated into the hot gas. 

Our set of initial condition s differs from that used in 
iBertone. Stoehr fc White] (|2005T ). who used the parametr i- 
sation for u wo and M w proposed bv lShu. Mo fc M ao (2005). 
In that paper, the wind mass loss rate and the wind initial 
velocity scale with the star formation rate as M w oc M^~ 0,71 
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Figure 1. K-band galaxy luminosity function at z = 0. Results 
from this work (sol id line) are comp ared to observ ational data 
from the 2dFGRS JColes et al.l l200ll) and 6dFGS lljones et al.l 
l2006h surveys (filled squares and circles respectively). The dashed 
line shows the corresponding luminosity function from the model 
of DLB07. 



and v wo oc °' 15 . The dependence of the wind initial ve- 
locity on the halo virial velocity in Eq. [8] can be expressed 
as v W o ^ 3.6i>„ ; ,. This is consi s tent w ith the parametri- 
sation of lOppenheimer fc Pavel |2006), wh o assume tha t 
Vwo ~ 3i>, r , r , in agreement wit h the findings of lMartinl l| 20051) . 
lOppenheimer fc Pavel l|2006l ) calculate the wind mass loss 
rate as M w oc v~}. This is equivalent to assuming q w = 1 
in Eq. [7] While our model assum es the same value of 
q w in Eqs. [7] and [8] the model of lOppenheimer fc Pavel 
(2006) requires two different values of q w . This implies 
th at the energy injection rate in winds Ei n in the model 
of lOppenheimer fc Pavel (|2006l ) scales with virial velocity 
as Ei n oc v vi[ , while the momentum input M w v vo is inde- 
pe ndent of v„ ir , as predicted for momen tum-driven winds 
bv lMurrav. Quataert fc Thompsonl (|2005l ). In our model we 
have assumed that the energy injection rate, and not the 
momentum injection rate, i s independent of halo proper ties, 
as assumed for example bv lSpringel fc Hernquistl (|2003h - 



3 THE ABUNDANCE OF GALAXIES 

In this Section we present results for the galaxy luminos- 
ity functions (LFs hereafter, Subsec. 13. the stellar mass 
function (Subsec. 13. 2p . and the morphology distribution of 
galaxies at z — (Subsec. 13. 3[) . We discuss galaxy colours in 
Subsec. 13.41 Throughout, we highlight the role of winds in 
suppressing star formation in faint galaxies, which is, per- 
haps, the most important achievement of our model. 



3.1 Galaxy luminosity functions 

Figs. [T] and [2] show the galaxy K-band and 6j-band LFs 
from our model (solid line) and the corresponding observa- 



Figure 2. As in Fig. [T] but for the 6j-band galaxy luminosity 
function at z = 0. 



tional data from t he 2 degree Field Galaxy Redshift Survey 
l|Coles et al.ll20qil, 2dFGRS herea fter) and the 6 degree Field 
Galaxy Survey ( Jones et ai]|2006l . 6dFGS hereafter) . In both 
figures, the dashed line shows the corresponding luminosity 
function from the model of ELB07. 

Our predictions differ from those of BLB07 both at the 
faint and at the bright end of the LFs. Our model predicts 
a shallower faint end slope than PLB07, bringing the abun- 
dance of the faint model galaxies closer to the observed val- 
ues. Although the agreement with observational data is still 
not perfect, it improves the results of PLB07 by a factor 
of two for galaxies fainter than Mk > —22. The ability to 
accurately reproduce the slope of the faint end of the galaxy 
LFs is peculiar to our model and is a consequence of our new 
feedback scheme. In addition, the slope of the galaxy LFs 
does not appear to be sensitive to the parameters regulating 
our feedback scheme. PLB07 predict steeper slopes for the 
faint end of the bj- and K-band LFs, which are not seen in 
observations. We have found that our model can reproduce 
these features only when feedback is particularly inefficient 
in dwarfs (i.e. q w = 0, or Vw and r\ have very low values). In 
this case, the energy injected in winds is not enough to over- 
come the gravitational potential of even the smallest haloes, 
leading to a severe overestimation of the number of faint 
galaxies. 

Figs. [T] and [2] also show that our model overpredicts 
the abundance of bright galaxies both in the bj- and in the 
K-band. This overprediction is a consequence of the short 
reincorporation time-scales obtained for relatively massive 
haloes. In our model, winds in haloes with M v ir > 10 12 
h~ Mq are unable to escape the gravitational potential and 
typically collapse after a few million years. The mass of the 
collapsing wind is instantaneously reincorporated into the 
hot gas and is again available for cooling, contributing to 
increase the luminosity of the central galaxy. This problem 
could partly be solved by assuming a higher effi ciency for the 
AGN feedback (20-30 per cent higher than in lCroton et al.l 
120061) or by increasing the reincorporation timescale of col- 
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Figure 3. Stellar mass function at z = 0. Numerical results for 
our model and f or D LB07 are compared to the observations of 
iBell et al.l ||2003|1 and lPanter et al.l J2004TI . 



lapsed winds in large haloes (as in |Pe Lucia et all |2004| ). 
However, since the aim of this work is to test the perfor- 
mance of a dynamical wind model in shaping the galaxy 
properties, assuming a non-instantaneous reincorporation 
time-scale defeats the main purpose of the model. 

3.2 Stellar mass function 

In this Subsection we focus on the stellar mass function of 
galaxies at z = 0, which we show in Fig. [3] Our re sults 
are compared to th e observations of IBell et al.l (|2003T ) and 
IPanter et all (|2004h and to the model of DLB07. We use the 
mass function results of IBell et all (|2003T ) obtained from the 
SDSS (/-band. Our model shows a better agreement with the 
observed stellar mass functions than DLB07, and in partic- 
ular it better reproduces the number density of faint galax- 
ies and of the brightest galaxies. However, both our model 
and the model of DLB07 underpredict the spatial density 
of M* galaxies at the knee of the distribution, for 10 10 ' 5 
ft -1 M < M* < 10 11 ' 5 ft _1 MQ. 

The stellar mass function should closely reflect the ri- 
band LF. Although at first sight the lack of galaxies around 
M* does not seem to appear in the K-band LF (Fig. (TJ, if 
we shifted the LF to lower luminosities by about 0.5 mag, 
we would recover exactly the same effect. This may be an 
indication that the galaxy luminosities are overestimated, 
or that, more likely, the galaxy formation model is unable 
to produce as many M* galaxies as there are in the real 
universe. The discrepancy between our predictions for the 
K-band LF and for the stellar mass function could also be 
partly due to a difference between the m ass-to-light ratios 
pre dicted by our model and those used bv lBell et al.l (|2003T j 
and lPanter et all l|2004l ). The high ag value used in the Mil- 
lennium simulation may further affect these results nega- 
tively. 

We believe that the underestimation of the abundance 
of M* galaxies is due to a combination of effects. In this 



Figure 4. Fraction of galaxy morphological types as a function of 
stellar mass. We define as ellipticals all galaxies with a dominant 
bulge component (Mb u i ge /M* > 0.7), as spirals all galaxies with 
< Mbuige/M* < 0.7, and as irregulars all galaxies with no bulge. 
M odel results are compared to the observational determination 
of lConselicel (|2006h . Thick coloured lines show the results of this 
work, while thin black lines are for the model of DLB07. 



range of stellar masses, the effects of galaxy and AGN feed- 
back overlap and as a result star formation may be sup- 
pressed more strongly than in other mass ranges. We find 
two possible ways to solve this problem. An overall increase 
in the efficiency of star formation or, alternatively, a decrease 
in the efficiency of feedback in galaxies as bright or brighter 
than M* helps to produce more galaxies at the knee, and 
in particular more early types. However, in the first case, 
much higher feedback is required to suppress the formation 
of dwarf galaxies and predict the correct slope of the faint 
end of the stellar mass function. On the other hand, reducing 
the efficiency of feedback in M * galaxies would imply that 
the efficiency to deposit supernova energy in winds depends 
critically on some unknown galaxy properties. In either case, 
there might be a connection between physical processes tak- 
ing place at scales not resolved by current simulations and 
the efficiency of star formation and feedback that we have 
not been able to include in our modelling of galaxy forma- 
tion. 



3.3 Galaxy morphologies 

In this section we analyse the morphological mix of our 
model galaxies at z = 0. Fig. [4] shows the fraction of galax- 
ies with different morphological types as a function of stel- 
lar mass. Data points in Fig. [4] show observational deter- 
minations from IConselicel (|2006t ). who analysed a sample of 
about 22000 galaxies with a vi sual morphological classifica - 
tion from the RC3 catalogue (|de Vaucouleurs et al.l [l99lh . 
As a proxy for the morphology of our model galaxies, we 
use the ratio between the bulge mass and the total stellar 
mass (r = Mbui ge /M*). In particular, we classify as ellipti- 
cals galaxies with more than 70 per cent of their stars in the 
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bulge (r > 0.7), as irregulars all galaxies without a bulge 
(r = 0), and as spirals all galaxies with < r ^ 0.7. The 
threshold value r = 0.7 is fixed arbitrarily, but results do 
not depend significantly on it. In particular, independently 
of the value of r, we find that ellipticals represent the dom- 
inant type at the highest stellar masses, while the smallest 
galaxies are predominantly irregulars. Most of the galaxies 
with 10 9 5 /i _1 M < M* < 10 11 h~ 1 M Q have a spiral mor- 
phology. A lower value of r would slightly decrease the frac- 
tion of spirals at intermediate stellar masses, while a higher 
value would slightly increase it. 

Fig.|3]shows that our model reproduces quite nicely the 
observational results. The main difference with the model 
predictions of DLB07 is found for the highest and lowest stel- 
lar mass galaxies. In the highest stellar mass bin our model 
does predict a few spiral galaxies, in agreement with obser- 
vations. The model of DLB07 does not predict any spiral in 
this mass range, even assuming a value of r = 0.6. This dif- 
ference is due to the recycling of gas in our feedback model, 
which is particularly efficient for galaxies in groups with 
Mvir > 10 13 /i _1 Mq. This results in galaxies with larger 
stellar masses to have larger disks in our model. We will 
discuss this feature of our model in more detail in Sec. [5] 
In the intermediate mass range, our model predicts a lower 
fraction of spirals than DLB07, in better agreement with 
observational results. In the lowest mass bins, our model 
predicts a non-zero fraction of spiral and elliptical galax- 
ies and a fraction of irregular galaxies lower than one. We 
note, however, that the morphological classification of model 
galaxies in this mass range might not be robust, as they re- 
side in haloes whose history can be followed back in time 
only for a few time-steps. 

In our model, the mass of the bulge is determined by 
two processes: mergers and dyn amical instabi li ties o f the 
disk. As expla i ned in detail in ICroton et alj |2006l ) and 
De Lucia et al.l (120061). th e disk instability is treated as in 
Mo. Mao fc White! (1998) and depends on the density of 
stars in the galactic disk. The outcome of a galaxy merger 
depends on the baryonic mass ratio of the merging galaxies. 
Our feedback scheme modifies the mass of gas and stars in 
galaxies (e.g. see Fig. [3| and so indirectly affects the mor- 
phological mix of model galaxies. 



3.4 Galaxy colours 

In this Subsection we briefly discuss the colours of galaxies 
and the colour-magnitude relationship. Fig.[5]shows the B — 
V colour of model galaxies as a function of stellar mass. The 
upper panel shows results from the model discussed in this 
paper while the lower panel shows resul ts from DLB07. Not e 
that this figure corresponds to Fig. 9 of lCroton et al. I (|2006l ). 
As in that work, the morphology of model galaxies is here 
defined on the basis of their bulge-to-total luminosity ratio 
in B-band. The vertical line indicate the limit above which 
the morphological classification can be considered robust. 

Fig. [5] indicates that there are several differences be- 
tween our model and the model of DLB07. While DLB07 
predicts a sharp colour bimodality, our model shows a much 
smoother distribution of galaxy colours, with few blue galax- 
ies and most of the galaxies residing on a red sequence 
which appears steeper and characterised by a larger scat- 
ter than the model results from DLB07. The suppression of 
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Figure 5. The B — V colours of model galaxies as a function of 
their stellar mass. The upper panel shows results from the model 
discussed in this paper while the lower panel shows result from 
DLB07. Galaxies are defined as "late types" and "early types" ac- 
cording to their bulge-to— total luminosity ratio in B-band. The 
vertical dashed line marks the limit where the morphological clas- 
sification of model galaxies can be considered robust. 



the colour-bimodality is a consequence of the wind model. 
We have verified that most combinations of parameters yield 
to a distribution similar to that shown in the upper panel 
of Fig. [5] A somewhat sharper bimodality can be obtained 
only in models with high values of the mass loss rate pa- 
rameter r). The results shown in Fig. [5] explain that our 
model predicts a lower (with respect to DLB07) abundance 
of red (B — V > 0.8), faint galaxies because the new feedback 
scheme leads to a steeper colour-magnitude relation. 
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Figure 6. Stellar metallicity as a function of stellar mass. The 
large dots and the shaded region represent the median and the 
16th and 84th percentiles of the distribut i on me asured for a sam- 
ple of SDSS galaxies by iGallazzi et~aH (|2005h . The solid and 
dashed lines show the same quantities for the model presented 
here and for the model of DLB07 respectively. 



4 CONSEQUENCES FOR THE CHEMICAL 
EVOLUTION OF GALAXIES 

In this Section we analyse how the stellar and gas metallici- 
ties of our model galaxies relate to their mass or luminosity, 
and compare our results to different observational measure- 
ments. The ability of a galaxy formation model to reproduce 
the observed correlations provides important constraints on 
how well the circulation of baryons in the different compo- 
nents is modelled. 

As explained in iKauffmannl (|l996h , the use of a surface 
density threshold for the star formation naturally reproduces 
the observed trend of the gas fraction as a function of galaxy 
luminosity. Our resul ts ar e very similar to th ose shown in 
|Pe Lucia et al.1 1|2004 ) and lCroton et alj (|2006f ) so we do not 
repeat them here. 



4.1 The metallicities of stars 

Fig.[6]shows our results for the metallicity of the stellar com- 
ponent of galaxies and compares them to the observational 
measurements by IGallazzi et al.l (|2005h and to the model 
results by DLB07. A ll model galaxies have been used here. 
IGallazzi et all <|2005h analyse about 44000 galaxies from the 
SDSS DR2 and find that the stellar metallicity increases 



with stellar mass up to M* 



10 1 



h 1 Mq. For larger- 



masses, the metallicity distribution flattens and saturates to 
a super-solar value. The stellar metallicity also flattens for 
Mi, < 10 9 ' 5 /i _1 Mq. Our model predicts a trend roughly 
similar to the observational data, although the distribu- 
tion of the simulated metallicities shows a smaller scatter. 
The model discussed in DLB07 exhibits an even smaller 
scatter and a shallower relation that appears to be in less 
good agreement with observational results. None of the two 



models exhibits the observed flattening for M+ < 10 
h~ 1 M(? l . Our predictions also broa dly agree with those of 
iMonaco. Fontanot fc Taffonil (|2007f ). 

The change in slope of the metallicity-mass relation in 
our model is a consequence of the feedback scheme. As we 
will show in Section [S] galaxies in haloes with M v j r < 10 12 
/i -1 Mq eject a larger fraction of their gas mass in winds and 
in the IGM than galaxies in larger haloes. This halo mass 
roughly corresponds to stellar masses of about M* ~ 10 10 ' 5 
h~ Mq. The high efficiency of mass ejection in small haloes 
strongly suppresses the formation of stars in these galaxies 
and therefore affects the amount of metals that can be pro- 
duced and locked into their stellar component, ultimately 
shaping the distribution of the stellar metallicities. We have 
experimented with different sets of parameters and verified 
that the flattening of the slope always appears where the 
efficiency of winds in ejecting gas from galaxies drops. 



4.2 Luminosity— gas metallicity relationships 

Fig. [7] shows the optical luminosity-gas metallicity rela- 
tionships in B-band and K-band and compar es the model 
result s wi th fits to observation a l data from Salzer et al.l 
l|2005l ) and lTremonti et alj (120041 '). iTremonti et al.1 (120041 ) es- 
timate the oxygen abundance in SDSS galaxies by fitting 
simultaneously the most prominent emission lines with a 
mode l for integrated gala xy spectra l|Charlot fc Longhettil 
l200lT ). ISalzer et al l (2005) use a sample of several hundred 
galaxies from the Kitt Peak Natio nal Observatory I nterna- 
tional Spectroscopic Survey (KISS, Sal zer et al.ll20001 ). They 
use two different metallicity calibrations: from Edmunds & 
Pagel (1984, hereafter EP) and from Kennicutt et al. (2003, 
hereafter KBG) . As can be seen from Fig. [7] different metal- 
licity calibrators give slightly different estimates of the oxy- 
gen abundance in nebular emission lines. Results from our 
model appear to be in v ery good agreement with the EP 
fit by ISalzer et~ai1 (|2005f l for Mb < -20 and M K < -24. 
At the brightest magnitudes, the luminosity-gas metallicity 
relation changes slope and model results appear to be closer 
to the fit corresponding to th e KBG calibrator. The esti- 
mate of ITremonti et aIT(|2004h intercepts the upper end of 
the distribution of the numerical results in B-band. 

To estimate the oxygen abundance of model galaxies, we 
assume a solar oxygen abundance of 12 + log 10 [0/H]q = 
8.76. The latest determinations of the solar oxygen abun- 
da nce range from 12 + log ln \Q / H]f^ = 8.66 as estimated 
by lAsplund. Grevesse fc Sauvall (J2006) using photospheric 



ab undances, to 12 + log in \Q/H]^ = 8.86 as estimated 



by iDelahave fc Pinsonneaultl (2006) on the basis of helio- 
seismology results. The higher value of 12 + log 10 [0/H]q 
wo uld bring our results in better agreement with the fit 
by ITremonti et all |2004), while the lower value would give 
a better agreement with the KBG fit. Independently of 
the normalization adopted, however, our predictions nicely 
match the observed slope of the luminosity-met allicity re- 
lationship in both B-band and K-band. Different combi- 
nations of model parameters shift the relation upwards or 
downwards by a few percent, but do not affect the slope or 
the scatter of the distribution. The dotted lines in Fig. [Jj 
show the corresponding relationships from the model of 
DLB07. When compared to our model results, this model 
predicts a shallower relation in the B-band and a non-linear 
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Figure 7. Gas metallicity as a function of the magnitude in 
B— band (upper panel) and in the K— band (lower panel). The 
black solid line shows the median of the distribution for model 
galaxies. The yellow marks the region including the 16th and 
84th percentiles of the di stribution. Observational determinations 
from [ Sa lzer et al] d2005T) (using two different calibrations) and 
iTremonti et al.l (1200417 are shown as dashed and dot-dashed lines 
respectively. 



relation in the K-band. Faint galaxies with A'Ib > — 18 and 
Mk > — 20 show the largest deviations from the observed 
metallicities. The model of DLB07 also features a smaller 
scatter than our model. 

Overall, the good agreement between the observed stel- 
lar and nebular metallicities and the predictions of our 
model is an indication that the dynamical feedback scheme 
for galactic winds is well integrated with the chemical evo- 
lution of galaxies. On the other hand, Figs. [6] and [7] indicate 
that the model of DLB07 may significantly overestimate the 
metallicities of faint galaxies, indicating perhaps an insuffi- 
cient ejection of metals in these galaxies. 



5 THE RECYCLING OF MASS AND METALS 
IN GALACTIC WINDS 

In this Section we focus our analysis on the recycling of 
metals in feedback processes. We trace the history of the 
deposition of metals into the IGM through cosmic time and 
investigate which haloes are the main sources of the metal 
enrichment of the IGM. Throughout, we highlight the main 
differences between our feedback scheme and that of DLB07. 

Feedback processes, and supernova feedback in particu- 
lar, determine the recycling of gas between episodes of star 
formation. Only a fraction of the galaxies in the simula- 
tion release enough energy in supernovae explosions to blow 
winds to distances comparable or larger than the virial ra- 
dius of the halo. As a consequence, not all winds can deposit 
mass and metals into the IGM. Most winds start blowing 
with such a little initial energy, compared to the energy re- 
quired to overcome the gravitational and pressure forces of 
the halo, that they collapse within a few time-steps. The role 
of these winds is important in regulating the star formation 
history of the central galaxies, because they determine the 
removal of the material reheated by supernovae explosions 
and its deposition into the hot halo gas. 

We assume that winds escape haloes when they reach 
the Hubble flow. When this happens, the wind mass and 
metals are deposited into the IGM and do not participate 
anymore in the evolution of galaxies in the halo. When winds 
reach the Hubble flow, the density field in which they ex- 
pand cannot be estimated analytically, unless the real 3- 
dimensional distribution of matter in the simulation is taken 
into account. We are therefore unable to trace the spatial 
distribution of these metals. In addition, our model neglects 
the possibility that metals mixed with the intergalactic gas 
can be incorporated at later epochs by nearby haloes or by 
the same parent halo accreting mass from a larger infall re- 
gion. 

Figs. [5] and [5] show how efficiently mass and metals, 
respectively, are ejected by winds and deposited into the 
IGM. The upper panel of Fig. [8] shows the average fraction 
of mass ejected by winds M out , defined as the sum of the 
mass currently in winds, plus the mass deposited into the 
IGM, Migm, over time. Similarly, the upper panel of Fig. 
shows the average fraction of the metal mass ejected by 
winds A/zout, again defined as the sum of the metal mass 
currently in winds, plus the metals deposited into the IGM, 
Mzigm, over time. Mbaryon is the total baryonic mass of 
the halo, while Mztot is the total metal mass in haloes and 
includes the metals in the gaseous phases and the metals 
locked in stars. The solid yellow lines are the fractions of 
ejected mass and metals in the DLB07 model at z = 0, 
respectively. The lower panel of Fig. [8] shows the average 
fraction of mass deposited into the IGM Migm, while the 
lower panel of Fig. [9] the average fraction of metal mass 
deposited into the IGM, Mzigm- Results are presented as a 
function of halo virial mass for z = 0, 1, 2, 3 and 4.5. 

The mass and metal ejection are most efficient in haloes 
with M v i r < 10 12 /i"'Mq. Winds from these haloes can eject 
up to 60 per cent of the baryons in the halo, of which more 
than half may have been deposited into the IGM. For haloes 
with 10 11 h~ 1 M P] < Af v ir < 10 12 h~ 'Mo , the DLB07 model 



predicts a mass ejection efficiency comparable to ours, but 
a significantly higher metal ejection efficiency. Larger haloes 
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Figure 8. Average fraction of ejected mass M ou t (upper panel) 
and of mass deposited into the IGM Mjgm (lower panel), as a 
function of halo virial mass. The ejected mass M out is the sum of 
the mass in winds, plus the mass deposited into the IGM Mjgm- 
The solid yellow line in the upper panel shows the fraction of 
ejected mass in the DLB07 model at z = 0. Results are shown at 
z = 0,1,2, 3 and 4.5. 



corresponding to groups and clusters can deposit at most 
a few percent of their mass and metals into the IGM, in 
agreement with the widely used assumption that the chem- 
ical evolution of clusters can be described at first order by 
a closed-box model. We will discuss the behaviour of the 
distribution of the ejected mass for M vir > 10 12 /i _1 M at 
the end of this Section. 

At z — 0, our model predicts a mass ejection efficiency 
significantly larger than that in the model used in DLB07 
for haloes with with M vir < 10 11 /i" 1 Mq. This difference is 
responsible for the lower abundance of dwarf galaxies pre- 
dicted by our model. The same haloes eject in our model 
on average twice the mass of metals ejected in the DLB07 



Figure 9. Average fraction of ejected metal mass Mz out (upper 
panel) and of metal mass deposited into the IGM Mugm (lower 
panel), as a function of halo virial mass. The ejected metals mass 
Mzout is the sum of the metals in winds, plus the metals deposited 
into the IGM Mzigm- Afztot is the total amount of metals in 
haloes. The solid yellow line shows the fraction of ejected metal 
mass in the DLB07 model at z = 0. Results are shown at z = 
0,1,2,3 and 4.5. 



model. This explains why our model produces on average 
dwarf galaxies with lower gas and stellar metallicities than 
the DLB07 model, as shown in Figs. [6] and [7] We note how- 
ever that haloes in this mass range contain typically less 
than a few hundreds particles and might therefore be af- 
fected by resolution. Results in this mass range should be 
interpreted with caution. The mass and metal ejection ef- 
ficiencies are only mildly dependent on redshift. However, 
metals are deposited into the IGM more efficiently with de- 
creasing redshift. 

Fig. [10] shows the cumulative distribution of the total 
metal mass density ejected by winds (upper panel) and of 
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Figure 10. Cumulative distribution of the ejected metal mass 
A^Zout (upper panel) and of the metal mass deposited into the 
IGM Mzigm (lower panel) per unit of comoving volume, as a 
function of halo virial mass. The solid yellow line in the upper 
panel shows the ejected mass in the DLB07 model at 2 = 0. 
Results are shown at z = 0, 1, 2, 3 and 4.5. 



not yet reached the Hu bble flow. Our findings are broadly 
consis tent with those of lBouche et alj (|2007l ). lBouche et al.l 
l|2007l ) investigated the metal budget in the IGM at z ~ 2.5 
and z ~ and found that galaxies with Lb < 0.3L B are 
responsible for about 90 per cent of the metal enrichment of 
the IGM. Our results for the mass density of ejected metals 
support their claim and we predict a similar amount of met- 
als ejected by galaxies per unit volume, as does the DLB07 
model. 



According to the upper panels of Figs. [5] and [5] haloes 
with Mvir > 10 12 ft _1 M can transfer up to 25% of their 
mass into a wind, before this is incorporated into the hot 
gas. This is an artifact of our feedback scheme. When a wind 
attempts to blow out of a massive halo, its initial velocity 
and the halo central gas density are high; this means that 
a much larger mass is accreted onto the wind from the halo 
gas than from the material outflowing from the star-forming 
galaxies. This effect becomes stronger with increasing halo 
mass, hence the trend of increasing ejection efficiency for 
Mvir > 10 12 /i _1 Mq. Since winds in massive haloes always 
collapse after a few timesteps, this feature simply describes 
the deposition of mass and metals into the hot cluster gas. 
The mass and metals ejected by haloes with M v i r > 10 12 
h~ Mq are therefore not relevant for the global mass and 
metal budget. This demonstrates to what extent our feed- 
back scheme is different from that of DLB07. This mass 
transfer do es not require an in termediate stage in the DLB07 
model (see lCroton et al.ll200rj for details). In DLB07, the hot 
gas is ejected only if the energy transferred to the hot gas by 
supernova explosions is larger than the increase of the ther- 
mal energy of the halo due to the deposition of the reheated 
mass. Since this never happens in haloes with M v i r > 10 12 
ft Mq, the hot gas is never removed from the halo. Hence 
the difference in the amount of gas being recycled in these 
haloes in the two models. The different treatment of feed- 
back in haloes with M v i r > 10 12 /i _1 Mq also explains why 
the model of DLB07 produces only red galaxies at the high- 
est stellar masses, while our model still predicts a few blue, 
late type galaxies (Fig. \%$. The continuous recycling of gas 
in groups with M vlr ~ 10 12 — 10 15 /i -1 Mq may favour the 
cooling of gas onto the central galaxies, which may accrete 
a sufficient amount of gas to form a disk. 



the metal mass density deposited into the IGM (lower panel) 
as a function of halo virial mass. Results are presented at 
z = 0, 1, 2, 3 and 4.5. The solid yellow line shows the results 
of DLB07 at z — 0. Mzigm and Mztot are calculated as in 
Fig.|9]and the mass density is expressed per unit of comoving 
volume. Fig. [TD] confirms the idea that the largest contribu- 
tion to the ejection of metals and to the metal enrichment 
of the IGM comes from galaxies in haloes with M v i r < 10 12 
ft Mq. At 2 = 0, the distribution of the ejected metals 
in the DLB07 model closely follows that predicted by our 
model for haloes with M vil < 10 12 /i _1 Mq. Fig. flOl also sug- 
gests that, while at z — most of the metals in the IGM are 
not gravitationally bound to haloes, at z = 2 — 3 most metals 
are still in the process of being ejected by winds and have 



To verify that our findings are genuine and do not criti- 
cally depend on the speficic values of the model parameters, 
we modify the power law index a w in eq. @ and ([8j and 
we consider models with q w = 0, 1 and 1.5 (for consistency, 
the value of the wind initial velocity Vw is also adjusted to 
keep the total wind energy roughly constant). The smaller 
the value of q w , the less mass is lost by smaller galaxies 
and the larger the mass injected in winds by larger galax- 
ies. However, we find that no matter the value of a w , haloes 
with Afvir < 10 12 /i _1 Mq are always the most efficient to 
deposit metals into the IGM, while larger haloes do not lose 
more than a few percent of their metals. This confirms that 
the dynamical feedback scheme we present here describes 
the evolution of galactic winds and their effects on galaxies, 
independently of the actual values of the model parameters. 
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6 CONCLUSIONS 

In this work, we have presented a new feedback scheme 
which is based on a dynamical treatment of the evolution 
of galactic winds in a cosmological context. The scheme has 
been incorporated in the Muni ch semi-analytic model for 
gala xy formation descri bed in |Pe Lucia fc Blaizot] (|2007r i 
and ICroton et alJ ([2006) and replaces previous implemen- 
tations of supernova feedback, but does not alter the AGN 
feedback model. In our scheme, the ejection and recycling 
of gas and metals are treated self-consistently and are gov- 
erned by the dynamical evolution of winds, rather than by 
empirical prescriptions. 

Our model contains a number of assumptions, as 
we have discussed above. These are: i) winds are 
described as spherically sym metric blastwaves, as in 
iBertone. Stoehr fc White] (120051 ). ii) The dynamics of winds 
is described as a two-stage process, e.g. a pressure-driven, 
adiabatic stage, followed by a momentum-driven snow- 
plough. An intermediate pressure-driven snowp lough stage 
is not included in our treatment of winds, iii) Winds sweep 
up all the gas they encounter along their path. This situa- 
tion correspond s to an entrainment fraction para meter e = 1 
in the model of lBertone. Stoehr fc White! (|2005| ). iv) Winds 
are linearly added together when haloes merge, without be- 
ing destroyed in the process, v) When winds collapse, the 
wind mass is instantaneously reincorporated into the hot 
halo gas. vi) We do not track the injection of energy into 
the hot halo gas when winds collapse and their mass is rein- 
corporated into the halo hot gas. vii) When winds reach the 
Hubble flow, their mass is permanently deposited into the 
IGM. viii) Winds are powered by the star formation taking 
place in all galaxies in the halo, including satellites, ix) The 
wind mass loss rate M w and the wind initial velocity 

^wo are 

described by Eqs. [7] and [8] respectively. 

Our results demonstrate that our dynamical scheme is 
a good alternative to empirical recipes for feedback and in 
particular it is able to reproduce the observed properties of 
dwarf galaxies with higher accuracy than previous models. 
Our main results can be summarized as follows: 

(i) The predictions for the the K-band and &j-band lu- 
minosity functions and for the stellar mass function are in 
good agreement with observations. In particular, the num- 
ber density of dwarfs is reproduced with greater accuracy 
than in the model presented in DLB07. 

(ii) The stellar mass-stellar metallicity and the 
luminosity-gas metallicity relationships are in good 
agreement with observations and significantly improve 
previous results. Again, the improvement is significant for 
dwarf galaxies. 

These results demonstrate that in our model the treat- 
ment of feedback in dwarf galaxies is significantly improved 
with respect to previous implementations. We argue that 
the main reason for this improvement is the possibility of 
estimating the amount of mass and metals that a halo can 
permanently lose and deposit into the IGM. This feature 
determines a stronger quenching of star formation in small 
haloes and therefore the suppression of the faint end of the 
stellar mass and luminosity functions of galaxies, and the 
shape of the mass-stellar metallicity and of the luminosity- 
gas metallicity relationships. 



Our model, however, also exhibits some weaknesses. 
Firstly, there is a tendency to overestimate the number of 
bright galaxies in the LFs, because of the continuous recy- 
cling of gas when winds recollapse in massive haloes. A pos- 
sible solution could be an increase in the efficiency of AGN 
fee dback (about 20-30 per cent more intense than suggested 
bv lCroton et~aill2006l ). but this would negatively affect the 
abundance of galaxies at the knee of the LFs. Our model 
does not solve the problem of the abundance of M* galaxies 
at z = 0, which continues to be underestimated (see also 
iBower et al.ll200"6h . Finally, perhaps the most severe short- 
coming of the model is that it does not correctly reproduce 
the colour distribution of galaxies. In particular, our model 
results do not exhibit the sharp colour bimodality observed 
for galaxies in the local Universe. This feature does not de- 
pend on the specific values of the model parameters and it 
can be alleviated only by assuming very inefficient feedback. 
However, models with inefficient feedback do not correctly 
reproduce the abundance of dwarf galaxies, nor their metal- 
licities. 

In the second part of our paper, we have investigated the 
efficiency of the mass and metal injection in winds and their 
deposition into the IGM. We have shown that in our model, 
galaxies in haloes with M v i r < 10 12 /i _1 Mq are responsible 
for most of the metal enrichment of the IGM, being the 
only objects that can blow winds out of their gravitational 
potentials. Larger haloes like groups and clusters lose only 
a small fraction of their mass and metals, and this occurs 
before most of the halo mass is accreted. 

Our feedback model contains the same number of free 
parameters ( i. e. 3) as the supernova feedback scheme of 
ICroton et all (|200tf ). once the depe ndence of the ejecte d 
mass on the halo virial velocity in ICroton et alJ (2006), 
which corresponds to our parameter a w , is taken into ac- 
count. The numerical calculation of the evolution of winds 
is only slightly more time-consuming than previous feedback 
schemes. However, it does require about 10-15 per cent ad- 
ditional memory to handle the wind variables and to store 
them on disk, if desired. 
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